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We studied the orbit of an electron revolving around an infinitely massive nucleus of a large 
classical Hydrogen atom subject to an AC electric field oscillating perpendicular to the electron’s 
circular orbit. Using perturbation theory in geometric algebra, we show that the equation of motion 
of the electron perpendicular to the unperturbed orbital plane satisfies a forced simple harmonic 
oscillator equation found in Lorentz dispersion law in Optics. We show that even though we did 
not introduce a damping term, the initial orbital position and velocity of the electron results to a 
solution whose absorbed energies are finite at the dominant resonant frequency ix = ujq\ the electron 
slowly increases its amplitude of oscillation until it becomes ionized. We computed the average power 
absorbed by the electron both at the perturbing frequency and at the electron’s orbital frequency. 
We graphed the trace of the angular momentum vector at different frequencies. We showed that at 
different perturbing frequencies, the angular momentum vector traces epicyclical patterns. 
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I. INTRODUCTION 

In standard optics texts, the position x of an electron 
of charge q and mass m under the time-varying electric 
field E = £e lut of light is given byPJSJ 

x+ Tx + coqX =-^—E, (1) 

2m 

where T is the damping coefficient, to o is the natural fre¬ 
quency of oscillation of the electron. The complex solu¬ 
tion x to this equation is shown by Akhmanov and Nikitin 
g] to be 


£ 

mix q — w 2 + iixY 


so that the power absorbed by the atom is 

„2 ,2p 

(P) = (qEx) = 


2m (wq — ix 2 ) 2 + <x 2 T 2 


( 2 ) 


|£| 2 , ( 3 ) 
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which in complex space becomes 

P = (£*x + £xe 2iut ) . (4) 

Notice that the damping term T makes the power ab¬ 
sorbed finite at the resonant frequency u = Xg. 

At present it is still not clear why an atom can be 
described as a simple harmonic oscillator subject to si¬ 
nusoidal electric field of light, e.g. what is responsible 
for the restoring force constant k = mix q and what is the 
cause of damping force Yx! 

In this paper we wish to show that a forced harmonic 
oscillator equation can arise for a large Hydrogen atom 
with a circular orbital of frequency wo subject to a lin¬ 
early polarized light of frequency ix whose corresponding 
wavelength 2'kc/ix is much larger than the electron’s or¬ 
bital radius r 0 , e.g. microwave frequencies, so that the 
phase of light is approximately the same at any point in 
the orbital path within a certain time period. That is, if 
the electron is initially in circular orbit in the xy- plane, 
the position s of the electron along the z-axis is given by 

s + u> 2 s = — — cos (cot) , (5) 

m 

as similarly given in Born and Wolf [5]. Even though 
this equation does not have a damping term, we shall 
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show that the energy absorbed at the resonant frequency 
u! = loq remains finite, provided we take into account the 
position r of the electron in 3D and use the vector form 
of the electrical energy dissipation expression pQ: 

P = c/E ■ f. (6) 

We shall show that (P) is finite at the resonant frequency 
uj = lo o even though there is no damping. 

In 1974, Bayfield and Koch experimentally studied the 
ionization of hydrogen atoms under microwave frequen¬ 
cies g5| . Since then, many tried to study the interaction of 
microwave radiation with classical hydrogen atom within 
the context of Rydberg atoms. [3 Some authors stud¬ 
ied the interaction with circularly polarized light [8Hl2]. 
while others such as Leopold na, Grosfeld and Friedland 
m and Neishtadt m focused on the linearly polarized 
case. 

For Leopold, his Hamiltonian is of the form: 

H{r,p) = ^p 2 - r _1 + zF max cosujt, (7) 

and the equations are solved using Monte-Carlo tech¬ 
niques. For Grosfeld and Friedland, their Hamiltonian 
is of the form: 

H = ip 2 — r _1 + Z^i cos’F, (8) 

where the frequency w(<) = d^z/dt and the authors used 
action angle variables. This Hamiltonian is the same one 
used by Neishtadt and Vasiliev, except that the latter 
authors used Delaunay elements. 

In our work, we shall not use the Hamiltonian ap¬ 
proach. Instead, we shall use the force equation 

fan 2 J* Q 

f =-r-p, - Ae 3 — E 0 cos(ut + ip), (9) 

to |r| J to 

and use linear perturbation theory to simplify the equa¬ 
tion to a simple harmonic oscillator equation in ([5]) for a 
motion perpendicular to the initial circular orbital plane 
of the electron. This method is simpler than those of the 
previous authors because the solution to the simple har¬ 
monic oscillator equation is well-known. Just as in the 
optical dispersion theory, we computed for the average 
power absorption by the atom and showed that it only 
depends on the ^-coordinate as in the standard theory. 

(P) T =- f Pdt. ( 10 ) 

T Jo 

We shall show that the orbit of the electrons at integral 
frequency ratios are similar to De Broglie waves, except 
that the oscillation is perpendicular to the electron’s or¬ 
bital plane. 

We shall divide the paper into six sections. Section 1 
is Introduction. In Section 2, we shall discuss the Ge¬ 
ometric Algebra formalism applied to planar rotations. 


In Section 3, we shall describe the unperturbed circu¬ 
lar orbit of the electron around the nucleus. After this, 
we shall introduce an oscillating electric field perturba¬ 
tion and derive the equations of motion of the electron’s 
oscillation perpendicular to its orbital plane, using the 
geometric algebra framework in our previous paper on 
Copernican epicyclical orbitsjTB]. In Section 4, we shall 
compute the electron’s orbital angular momentum and 
determine its limiting form at the resonant frequency. In 
Section 5, we shall compute the electrical power dissi¬ 
pation of the electron and plot the results for different 
values of the ratio between the orbital and light frequen¬ 
cies. We shall show that the average power, either over 
the perturbing or orbital period, is approximately simi¬ 
lar to the standard absorption resonance curve with finite 
peak. Finally, we graph the angular momentum of the 
electron at different frequency ratios and show that the 
angular momentum vector traces epicyclical patterns. 

II. GEOMETRIC ALGEBRA 
A. Scalars, Vectors, Bivectors, and Trivectors 

In Clifford (Geometric) Algebra Clz,a, also known as 
the Pauli Algebra, the product of the three unit vectors 
e!,e 2 , and e 3 satisfies the orthonormality relation P1HIS] 

&j^k 4“ &k&j — 2(5 jk, (11) 

where 5jk is the Kronecker delta function. In other words, 
the square of the length of the vectors is equal to one and 
the product of two perpendicular vectors anticommute. 

Let a and b be two vectors spanned by e 2 , and 
e 3 . We can show that their product satisfies the Pauli 
identity[20] 2b 

ab = a • b + i{ a x b), (12) 

where i = e!e 2 e 3 is the unit trivector which behaves like 
an imaginary scalar that transforms vectors to bivectors. 
The Pauli identity states that the geometric product of 
two vectors is equal to the sum of their scalar dot product 
and their imaginary cross product. 

B. Exponential Function and Rotations 

Let ie 3 d be the product of a bivector ie 3 = eie 2 with 
the scalar 6. Since the square of ie 3 0 is negative, then 
the exponential of ie 3 6 is given by Euler’s theorem 

— cos 0 -f- ze 3 sin 9. (13) 

From this we can see that 

cos8= ^(e te3d + e~ le30 ), (14a) 

sin0 = —-—(e* 63 ® — e - *® 3 ®), (14b) 

2 *e 3 
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FIG. 1: Rotation of ei and e 2 about e 3 counterclockwise by 
an angle 9 


which are the known exponential definitions of cosine and 
sine functions. 

Multiplying Eq. (131 by ei, e 2 , and e 3 , we obtain 


eie*® 3 ® = ei cos 0 + e 2 sin0 = e -M33e ei, (15a) 

e 2 e* e3 ® = e 2 cos 6 — e! sin 6 = e~ te30 e 2 , (15b) 

e 3 e lB30 _ e3 cog q _|_ e3 j e3 s i n Q = e le:>e e 3 . (15c) 


Notice that eie*® 30 is a rotation of ei counterclockwise 
about e 3 by an angle 8 , while e 2 e !e3S is a rotation of e 2 
counterclockwise about the same direction and the same 
angle. Notice, too, that the argument of the exponen¬ 
tial changes sign when e 3 or e 2 trades places with the 
exponential, while e 3 commutes with the exponential. 

A vector a in 2D can be expressed in both rectangular 
and polar forms: 


In general, we have the following relations: 


eia = a*ei, (20a) 

e 2 a = a*e 2 , (20b) 

e 3 a = ae 3 . (20c) 

That is, e 3 and e 2 both changes the complex number a 
to its conjugate a* after commutation, while e 3 simply 
commutes with a mrnrn- 


III. LIGHT-ATOM INTERACTION 


A. Unperturbed Electron Orbit 


Classically, the position r of an electron of mass m 
and charge — q as it revolves around a massive proton of 
charge q is given by Coulomb’s law: 


kq 2 r 
m |r| 3 ’ 


( 21 ) 


where k is the electrostatic force constant. We claim that 


a solution to Eq. (21) is given by 


ro = eif 0 V>o, 


( 22 ) 


where 


f 0 = r 0 e ie3Vo (23a) 

= e ie ^ot (23b) 


are the complex amplitude and the rotation operator, 
respectively. Substituting these back to Eq. (22), we get 

r 0 = e 1 r 0 e ie3( " J ° t+Vo) , (24) 


which yields 


a = a x e 3 + a y e 2 = ae 3 e 


ie 3 9 


(16) 


ro = ep'o cos (u 0 t + <p 0 ) + e 2 r 0 sin (w 0 t + <p 0 ), (25) 


Expanding the exponential using Eq. (15a) and separat 


ing the e 3 and e 2 components, we arrive at the standard 
transformation equations for polar to rectangular coor¬ 
dinates: 


x = a cos 9, 
y = asm. 9. 


(17a) 

(17b) 


We may also factor out e 3 in Eq. (16) either to the left 
or to the right to get 


a = e 3 d = e 3 (x + ie 3 y ) = e 


ie 3 6 


i ae 


a = a*ei = (x — ie 3 y)e 3 = ae le3t ei. 


(18a) 

(18b) 


Factoring out e 3 yields the definition of the complex num¬ 
ber a and that of its complex conjugate a*: 


a = a x + ie 3 a y = ae*® 3 ®, 
a* = a x — ie 3 a y = ae~ le30 . 


(19a) 

(19b) 


after expanding the exponential and distributing e 3 . 
Equation (22) states that the electron moving around the 


proton in circular orbit of radius r 3 with angular velocity 
wo and rotational phase angle <po. 

To verify that Eq. (22) is indeed a solution to the 
Coulomb’s law in Eq. (21), we first compute the first 


and second time derivatives of Eq. (221: 


r = eif 0 ze 3 w 0 V>o, 
r = -eifoWQ^o. 


(26a) 

(26b) 


Now, substituting Eqs. (22) and ( |26b ) to the Coulomb’s 
law in Eq. (21), we obtain 


OG 2 = 


kq 2 


„3 5 


(27) 


after cancelling out e^o^o. Equation (27) is the familiar 
circular orbit condition. 
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FIG. 2: Uniform circular motion of an electron with a distance 
r o from the nucleus. The orbital angular frequency is to and 
the phase angle is i p. 


F P 



FIG. 3: The Coulomb force F c and the perturbing force F p 
on an electron moving with velocity vq and radius r o 


B. Perturbation by an Oscillating Field 


Suppose that the electron is subject not only to the 
Coulomb force — qE c due to the proton, but also to the 
force — qEi p due to an oscillating perturbing field. The 
equation of motion of the electron then becomes 


rnr = — qE c — AgEp, (28) 


where A is a perturbation parameter that shall later be 
set equal to unity. More specifically, we write 


r 


kq 2 
m 



—E 0 cos {tot + ip), 
m 


(29) 


where Eq,w, and tp are the amplitude, angular frequency, 
and phase of the perturbing electric field. Our aim is 
to determine the position r of the electron that satisfies 
Eq. ((29]). 

To find the solution to the perturbed orbit equation 
in Eq. (29), we assume that the solution is a sum of the 
electron’s unperturbed circular orbit in Eq. (22) and a 
slight perturbation s perpendicular to this orbit. So we 
write 


The first and second time derivatives of r are 


r = r 0 + Ari = + Ae 3 s, (31a) 

? = ? 0 + A? i = -eifow^o + Ae 3 s. (31b) 


Equation (31b I shall take care of the left side of Eq. 

To expand the right-hand side, we need first to take 
the square of the position vector r in Eq. (30) and retain 
only the terms up to first order in A: 


r 2 = r„ + 2A(r 0 ■ iq). (32) 


Since r 0 lies on the unperturbed orbital plane of the elec¬ 
tron xy plane and iq = e 3 is perpendicular to this plane, 
then r 0 • iq = 0, so that Eq. (32) reduces to 


r 2 = rl = e 1 fo'i/’o e^o^o = e 1 e 1 f2V>2n)V’o 
= rgro = rl, 


(33) 


where we used the definitions of r o and ip q in Eqs. (23a) 
and (23b). Thus, |r| = ro, so that 


■j—rq = - 3 eifoV’o- 

| r | 3 r 3 


(34) 


Equation (34) shall take care of the Coulomb term on the 
right side of Eq. (|29|) . 


Now, substituting Eqs. (31b) and (34) back to equation 
of motion in Eq. (291, we obtain 


-eif 0 w 2 ^o + A e 3 s = -to^eif^o + Ae 3 s) 


— Ae 3 — Eq costut + tp), (35) 

TO 


where we used the circular orbit condition in Eq. (27). 
The term zeroth order in A cancels out, so we are left 
with the term first order in A. Hence, 


s + WqS = — — Eq cosfcot + tp), (36) 

TO 

after rearranging the terms. Notice that Eq. ( |36| ) is a 
simple harmonic oscillator equation with sinusoidal forc¬ 
ing, which is the standard model for classical light-atom 
interaction. 


C. Solving the Forced Harmonic Oscillator 
Equation 


Let §h and s p be the homogeneous and particular so¬ 
lutions of the Eq. (36). That is, 


s — Sh + s p 


(37) 


and 


Sh + u]Q§h = 0, (38) 

Sp + u;gSp =- Eq cos(wf + tp). (39) 


r = r 0 + Ari = eifoV’o + Ae 3 s. 


(30) 
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The solution to the homogenous equation in Eq. (38) 
is a sum of a sines and cosines: 


Sh = Chi cos(w 0 t) + c h 2 sin(w 0 t), 


(40) 


where c^i and Ch 2 are scalar constants that will be de¬ 
termined from the boundary conditions. On the other 
hand, the solution to the particular equation in Eq. (39) 
is of the same form as the perturbing held: 


s P = c p cos(wi + p), 


(41) 


where c v is a scalar constant. Substituting Eq. (41) back 
to the particular equation in Eq. (39) and solving for c p , 
we get 


qE 0 


1 


mw§ (a 2 — 1) ’ 


where 


w 


a = 


UJ 0 


(42) 


(43) 


is the ratio of the perturbing frequency w to the electron’s 
orbital frequency ojq. Hence, 


1 


Sr, — 


m«g (a 2 — 1) 


Eq cos(awot + p)- 


(44) 


Adding the homogenous solution Sh in Eq. (40) to the 
particular solution s p in Eq. (41) yields the total solution: 


s = c h i cos(w 0 t) + c h2 sin(w 0 t) 

q 1 


+ 


mojQ (a 2 — 1 ) 


Eq cos(auiot + p). (45) 


Its time derivative is 

s = -ChiLOo sin(w 0 t) + c h2 u)Q cos(w 0 t) 

n rx 

-Eq sin(awot + p)- 


(46) 


mwo (a 2 — 1) 

To determine the unknown constants Chi and Ch 2 , we 


hrst substitute the expressions for s and s in Eqs. (45) 


and (46) back to the expressions for the position r and 


velocity r in Eqs. (30) and (31a) to get 


r = eifo^o + e 3 (chi cos(w 0 t) + c h2 sin(w 0 t)) 


+ e 3 “ 2/2 

mr, it i 


1 


(a 2 - 1) 


Eq cos(awot + p ), 


r = eif 0 ie 3 w 0 ^o + e 3 (-c hl u} 0 sin(w 0 t) + c h2 u o cos(w 0 t)) 

e 3 ~ 7 2 ° ' E 0 sm(au;ot + p), (47b) 

tow o {a z — 1) 

after setting the perturbation parameter A = 1. If we 
assume that at t = 0, the electron is in its unperturbed 
circular orbit around the nucleus, then 


r(0) = eir 0 , 
r(0) = eir 0 *e 3 w 0 . 


(48a) 

(48b) 


Substituting these to Eqs. (47a) and (47b), and setting 


t = 0, we arrive at the expressions for the parameters Chi 
and c h2 : 


1 


Chl = - 


Ch2 = 


toWq (a 2 — 1) 


Eq COS ip, 


towq (a 2 — 1) 


Eq sin tp. 


(49a) 

(49b) 


Substituting Eqs. (49a) and (|49b ) back to the expres¬ 
sion for the position r in Eq. (47a), we get 


r = eir 0 ^o + e 3 


qE 0 


'toWq (a 2 — 1) 

+ sin ip sin(Wot) + cos(awot + tp)). 


(— cos^cos(wot) 


(50) 


Using the identity for the cosine of a sum of two angles, 
Eq. (l50|) reduces to 


r = eifo^o + e 3 


qE 0 


1 


tow § (a 2 — 1) 


(cos (aw 0 t + ip) — cos(w 0 t + p))- 
Its time derivative is 

i . qE 0 l 

r = eir 0 /e 3 w 0 ^o + e 3 -- —r x 

towq \or — 1 ) 

(—a sin(awot + ip) + sin(wo t + ip), 


(51) 


(52) 


where we used the definition of a = w/w q. Equations (51) 


and (52) are the position and velocity of the electron 
initially orbiting at radius r 3 , angular frequency wo, and 
phase ipo, and perturbed by an oscillating electric field 
with amplitude Eq, frequency w, and phase ip. 


D. Orbit Equations and Limiting Conditions 


To convert Eqs. (51) and (52) into rectangular coor¬ 


dinates, we use the expansions in Eq. (15a) and (15b), 


together with the identity eiie 3 = e 2 to arrive at 

x = ro cos(wot + tpo), (53a) 

V = tq sin(w 0 t + <po)t (53b) 

qE 0 1 


z = 


tow 2 (a 2 — 1) 


(cos(awof + tp) — cos(wot + p)), 

(53c) 


(47a) and 


x = -r 0 w 0 sin(w 0 f + p 0 ), 
y = r 0 w 0 cos(w 0 t + ^o), 


(54a) 

(54b) 


z = *1— - -— (—asin(aw 0 t + p) + sin(w 0 t + p)). 

mu>o ( a— 1) 

(54c) 


Equations (53a) to (53c) are the equations for plotting 


the orbit of the electron as a function of time. Equa¬ 
tions (54a) to (54c) are for plotting the corresponding 


velocities. 
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— unperturbed orbit 


— unperturbed orbit 

- o=0 


- q=0.5 




FIG. 4: The unperturbed orbit lies flat along the ei — e 2 plane. When a = 0, or when the electric field is constant in time, the 
orbit slants. When a —» oo, more waves are observed. 


When the perturbing frequency u> = 0, corr espon ding 
to a = 0, the expressions for z and i in Eqs. (53c I and 
(54c) reduces to 


2 

z 


qEp 

mojQ 

qEp 

mojo 


(cos If — COs(bJpt + ¥>)), 
sin(u;ot + <f)- 


(55a) 

(55b) 


If ip = 0, the perturbing field is E = E Q e 3 , so that 


z = — ^ (1 — cos(u>ot)), (56a) 

mu>o 

z = — sin(woO- (56b) 

mwo 

On the other hand, if ip = tt, the perturbing field is E = 
E 0 e 3, so that 


z = g °o (1 + cost (57a) 
mu> o 

z = sin(wot)- (57b) 

mu> o 

These are the behavior of the electron’s orbit along the 
direction when the perturbing electric field is constant, 
also known as the DC electric field. 

Now, when the perturbing frequency ui = ujp, corre¬ 
sponding to a = 1, the field resonates with the electron’s 


orbit. The only terms affected are z and z in Eqs. (53c) 


and (54c I. Since both their numerators and denomina¬ 
tors approach zero as a —>• 1, we apply L’hopital’s rule by 
differentiating the numerators and denominators prior to 
evaluation of the limits: 

lim z = lim , ( 58a ) 

ct->l TOWq a-> 1 y 2 a J 


lim i = lim 

a->l mwo a-*- 1 


qE 0 


lim 


sin(awoi + f) 
2a 

f aojpt cos(ojt - 




mwo “-* 1 \ 


2a 


Hence, 


lim x = 

a —>-1 

lim z = 

ot—y 1 


qEp 

2mujp 

qEp 

2mujp 


t sin(auJot), 


(58b) 


(59a) 


(sin(wot + ip) + uipt cos(ojpt + if)). 


(59b) 


Notice that the amplitude of the oscillations along z and 
its corresponding velocity are linearly increasing in time. 
Once the amplitudes of the oscillations becomes so large, 
our perturbation approximations breaks down. Thus, our 
theory cannot really say what happens during ionization 
or whether ionization will really happen at all at the res¬ 
onant frequency. (See Figs. [4] and [5]) 
































7 



FIG. 5: Height z of the electron from its unperturbed circular 
orbit with respect to time. 


IV. ANGULAR MOMENTUM 
A. Product Form 

Let us compute the product of t he position r in 
Eq. (30) and its velocity r in Eq. (|31a[) : 


rr = (eif 0 t/>o + Ae 3 s)(eiroze 3 woV’o + Ae 3 s). 
Distributing the terms, we get 
rr = e 1 e 1 fl^r 0 ie 3 uj 0 'ip 0 

+ e 3 ei(sf 0 fe 3 u;o^o - sr^) + e 3e 3 ss, 


(60) 


(61) 


after setting the perturbation parameter A = 1. Since 
i = eie 2 e 3 and ie 3 = eie 2 , then Eq. (61) reduces to 


rr = ie 3 iu 0 rQ — ie 1 sfg(UoV’g — ie 2 sfg^>g + ss. (62) 


Separating the scalar and bivector parts of Eq. (62), we 
arrive at 


r • r = ss, (63a) 

„2 „(63b) 


r x r = e 3 w 0 ro - eisf 0 w 0 ^o - e 2 sr 


after factoring out the trivector i in the second equation. 
Multiplying Eq. (63b I by the electron’s mass m yields 


the the electron’s angular momentum: 

L = mr x r = e 3 w 0 ro - eisfowo^o - e 2 sfg^o, 


(64) 


where s and s are the e 3 components of the electron’s 
position r and velocity r in Eqs. (51) and (52): 


qE 0 


1 


?nujg (a 2 — 1) 


qE 0 


1 


muiQ ( a 2 — 1 ) 


(cos(awot + — cos(uiot + ip)), 

(65a) 

(—a sin(awgt + ip) + sin(u;gi + p ). 

(65b) 


Using the definitions fg = rge lG3¥3 and tpo = e 




Eq. (64), and separating the ei, e 2 , and e 3 components, 


we arrive at 

L 1 = —mwosro cos(wot + tpo) — msr g sin(wot + <po), 

(66a) 

L 2 = — mijJosro sin(wot + tpo) ~ msro cos(wot + ipo), 

(66b) 

L 3 = mw 0 »o, (66c) 

which are the parametric expressions for the angular mo¬ 
mentum in rectangular coordinates. 

B. Harmonic Form 

The vertical oscillation s and its derivative s in 


Eqs. (65a) and (65b) may be expressed in exponential 
forms: 


qE 0 


1 


O _ _ ( p ie 3 (au 0 t+ip) I -ie 3 (aui 0 t+ip) 

2 toWq (a 2 - 1) 1 + 

_ e *e 3 (w 0 i+v) _ e -ie 3 (uJot+ip) j (67a) 

— ie 3 (aui 0 t+ip) 


g — g-S'O ^ e 3 ^ ac ie 3 (au>ot+ip) 


2mwo (a 2 — 1) 
These may be rewritten as 


e ie 3 (u] 0 t+<p) _ e 


ae 

-ie 3 (u 0 t+ip ) \ 


(67b) 


—a. 


1 E 0 1 

- #o ~ V > 0 ’ 1 )> 

a = ( J ° , n 3 ,, (—afj-ip o + “ 


2muo (a 2 — 1) 


+ #o - V*i >o 1 )> 


where 


fj = e ie3V> . 


(68a) 

(68b) 

(69) 


noting that e 2 ie 3 = — ei, we obtain 


Substituting Eqs. (|68aJ) and ( |68b| back to Eq. (|64) and 

(70) 


L = e 3 mojor 3 — e 


qE 0 


1 


2w 0 (a 2 - 1) 


i>L, 


where 


ipL = f/roTpo + V*roipo ~ ~ fj*ro 


o 1 + af)*fg^ 0 a 1 +Wo-Tj*foip 0 ' 2 . 

(71) 

Expanding the terms ofipL into exponential form, we get 

^j) L _ r ^ e ie 3 {.{a+l)u] 0 t+ip+Lfi 0 ) e ie 3 ((-Oi+l)ui 0 t-y+<p 0 ) 

— e ie 3{^at+v+Vo) _ e ie 3 (-(p+<p 0 ) 

_ ae ie 3 {(a-l)u> 0 t+<p-tp 0 ) + ae ie 3 (-(a+l)u 0 t-ip-v 0 ) 


e ie 3 (ip-ip 0 ) _ e ie 3 (- 2 ui 0 t-<p-tpo) 


), 


(72) 































Substituting the result back to Eq. (70) and separating 


the ei, e 2 , and e 3 components, we arrive at 
qE 0 r 0 1 


L i — — 
L2 = — 


2 w 0 (a 2 — 1 ) 
qE 0 r 0 1 
2 w 0 (a 2 - 1 ) 
L 3 = mwoTo, 


7i, 

(73a) 

72, 

(73b) 


(73c) 


where 


71 = (1 + a) cos((a + l)w 0 t + + <A)) 

+ (1 — a) cos((a — l)w 0 t + <p — ipo) 

— 2 cos (2 io 0 t + p + po), (74a) 

72 = (1 - a) sin((a + l)w 0 t + p + p 0 ) 

— (1 + a) sin((a — l)w 0 1 + p — po) 

+ 2sin(v? — tpo). (74b) 

Thus, since a = w/w 0 , we see that the orbit of the tip of 
the angular momentum vector L is a linear combination 
of circular motions with the following orbital frequencies: 

WL = {±(w + w 0 ),±(w - w 0 ), ±2w o ,0}. (75) 


V. POWER AND ENERGY ABSORPTION 
A. Work-Energy Theorem 

The work-energy theorem states that 


' , 1 2 1 2 
t ■ dr = -mv — -mvg . 


This may be rewritten as 

rt 
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1 2 1 2 

Pdt= -mv 2 - 2 mv o ’ 


where the power P is defined as 

P = F • r. 


(78) 


(79) 


(80) 


That is, the integral of the power P expended by a force 
F acting to move a mass m from time t = 0 to t is equal 
to the change in the mass’s kinetic energy between these 
times. 

In our model, there are two forces acting on the elec¬ 
tron: the Coulomb force F c and the perturbing force F p . 
But since the sum of these two forces is mr, as given in 
Eq. (f2Sl) , then the left side of Eq. (l79|) becomes 


C. Limiting Conditions 

In the DC field limit, a —» 0, so that 


P = mr • r. (81) 

We shall use this equation to compute the power ab¬ 
sorbed by the atom. 


qE 0 r 0 

= -X-(cos (w 0 f + p + p 0 ) 

ZUJo 

+ COs(— w 0 t + p — Po ) 

— 2cos(2wot + p + vpo))> (76a) 
qE 0 r 0 . 

= -x-(sm(w 0 t + p + po) 

2 w 0 

- sin(—w 0 f + p — po) 

+ 2sin(</j — po))- (76b) 

On the other hand, in the resonance frequency limit, a -A 
1 , both the numerators and denominators of L\ and L 2 
approach zero, so that we apply the L’hopital’s rule: 

lim Li = — 1 ° ° (cos( 2 wgt + ip + po) 
a ->0 4wo 

- 2w 0 f sin(2w 0 t + <p + po) - cos(y> — Po)), 

(77a) 

lim L 2 = — j , ° ° (— sin( 2 w 0 t + p + po ) 
a —>0 4Wo 

— sin(y> — ipo) — 2wo t cos(<p — po))- (77b) 

Notice that at the resonant frequency w = w 0 , the L\ 
and L 2 components of the angular momentum increases 
in time; in our perturbative approximation, the atom will 
be ionized. 


lim L\ 

CL —^0 


lim L 2 

ot —^0 


B. Power: Product Form 


To evaluate the left side of Eq. (81), we firs t mu ltiply 
the expressions for r and r in Eqs. (31b) and (31a): 


rr — (—ejygWQ^o + Ae 3 s)(e 1 fo*e 3 Wo^ , o + Ae 3 i) (82) 
Distributing the terms, we get 


rr = —e 1 e 1 fgWg'!/jQf 0 je 3 wo'i/'o + e 3 eisfoie 3 w 0 ^g 

,2 


eie 3 f 0 w 0 V’ 0 'S + e 3 e 3 ss, 


(83) 


after setting A = 1. Since i = eie 2 e 3 and ie 3 = eie 2 , 
then Eq. (83) reduces to 


rr = — ie 3 rgWg — ieiif oWg^g + *e 2 srgWg^o + si. (84) 
Separating the scalar and bivector parts, we get 

r • r = si, (85a) 

r x r = -e 3 rgWg - eiifoWo^o + e 2 ifoWgV)o, (85b) 

after factori ng ou t i in the second equation. 

Equation (85a) leads to a very simple expression for 


the power absorbed by the atom: 

P = mi ■ r = mss. 


( 86 ) 
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FIG. 6: Projection of the tip of the angular momentum vector L in the xy plane for different values of a. 


Taking the time derivative of s in Eq. (65b), 


s = -—-tv (—a 2 cos(awo t + </?) + cos(u>ot + </?)), 


m (a 2 — 1) 


(87) 


and substituting this and that of s to Eq. ( 86 ), we get 


P = 


1 ;; 

2 mwj (a 2 — l ) 2 

(—a 2 cos(awot + ip) + cos(wo t + tp)) x 
(—a sin(au;of + ip) + sin(wot + <p)). 


( 88 ) 


We can show that this is equivalent to 
~ 2^2 2 

P = - - - --vr (a 3 sin(2aw 0 t + 2 p) 

zmujQ (or — 1 y v 

— (a 2 + a) sin((a + l)wof + 2 (p) 

+ (a 2 — a) sin((a — l)wot) 

+ +sin(2wot + 2ip)), (89) 



FIG. 7: This is the graph of the average power over the per¬ 
turbing period 


which is the desired harmonic form of the power P ab¬ 
sorbed by the atom. Notice that power is not constant 
but fluctuating in time. 
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C. Average Power over Perturbing Period 


Let us define the average power over the perturbing 
period as 


where 



Now, let us take the time average of the power P in 
Eq. (98) over the perturbing period: 


< P) ' = (STTip PM W + 

-(a 2 + a)(sin((a + 1 )w 0 f + 2 <p)) T 
+ (a 2 - a)(sin((a - l)w 0 t))r 
+(sin(2a;ot + 2<p)) T ), (92) 



a 


FIG. 8: This is the graph of the average power over the orbital 
period 


where 


D. Average Power over Orbital Period 


(sin(2awot + 2 (p)) T 
= 0 , 

(sin((a + 1 )uj 0 t + 2 ip)) 

1 a 


(93a) 

(cos((l + l/a)2-7r + 2 ip) — cos(2 ip)) 

(93b) 

(cos((l — l/a)27r) — 1), (93c) 


2-7T (a + 1) 

(sin((a - l)w 0 t)) T 
1 a 
2ir (a — 1) 

(sin(2wot + 2( p)) T 

= -!-a(cos(47r/a + 2tp) — cos(2<p)). 

Z 7 T 


(93d) 


Substituting these back to Eq. (98), we get 

J2 rp2 

(P)r= 9 ° 


1 


47rmwo (a 2 — l) 2 
(—a 2 (cos((l + l/a)27r + 2ip) — cos(2(/j)) 
— a 2 (cos((l — l/a)27r) — 1) 

+ a(cos(47r/a + 2 ip) — cos(2 ip)). 

If tp = 0, we have 

fEt 1 .. 


(94) 


(P)r = 


47rjncuo (a 2 — l) 2 
(—a 2 (cos((l + l/a)27r) — 1) 
— a 2 (cos((l — l/a)27r) — 1) 
+ a(cos(47r/a) — 1)). 


(95) 


Equation (951 is graphed in Fig. [7j Notice that despite 
the small oscillations in the interval a = 0 and a = 1, 
the power absorption curve is similar to that in Lorentz 
dispersion theory. 


Let us define the average power over the perturbing 
period as 


1 f T ° 

(P)t o = — Pdt , 

To Jo 


where 


to = 


2 tt 
u o' 


(96) 


(97) 


Now, let us take the time average of the power P in 
Eq. (f98l) over the orbital period: 


„2 Z7>: 

(P)r, = “ E ° 


r (a 3 (sin(2awot + 2 ip)) 


_ _ 1 

T ° 2mwo (a 2 — l) 2 ' 

- (a 2 + a)(sin((a + 1 )w 0 t + 2 ip)) To 
+ (a 2 - a)(sin((a - l)ui 0 t)) To 

+ (sin(2w 0 t + 2ip)) To ), 


TO 


(98) 


where 


(sin(2awot + 2 ip)) To 

1 cos(47ra + 2 ip) — cos(2 ip) 
~2n 2a ’ 

(sin((a + l)w 0 t + 2 ip)) To 

1 cos(27r(a + 1)) — cos(2(/j) 
2-7T a + 1 

(sin((a- l)w 0 t)) To 

1 cos(27ra) — 1 
27r (a — 1) 

(sin(2wot + 2tp)) To 
= 0 . 


(99a) 

(99b) 

(99c) 

(99d) 
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Substituting these back to Eq. (98), we get 

1 


q L ( 


(P)r o = 


T ° 8Trmui 0 (a 1 2 3 4 5 6 7 8 * 10 11 — l) 2 

(—a 2 (cos(47ra + 2ip) — cos(2 ip)) 

— 2a(cos(2ir(a + 1)) — cos(2(p)) 

- 2a(cos(27ra) - 1)). (100) 


If ip = 0, we have 

q 2 E 2 -1 


r (a 2 (cos(47ra) — 1) 


^ ^ T ° 87rma;o (a 2 — l) 2 ' 

— 2a(cos(27r(a + 1)) — 1) — 2a(cos(27ra) — 1)). 

( 101 ) 

Equation ( |101 ) is graphed in Fig. [8] Notice that unlike 
in Fig. [TJ there are periodic oscillations after a = 2. 


VI. CONCLUSION AND RECOMMENDATION 

We modelled the classical Hydrogen atom as an elec¬ 
tron revolving in circular orbit around an immovable pro¬ 
ton subject to the Coulomb force. We subjected this 
atom to a perturbing oscillating electric field perpendicu¬ 
lar to the electron’s initial orbital plane. We showed that 
the resulting equations of motion of the electron along 
the axis of the perturbing electric field is similar to that 
of a simple harmonic oscillator with sinusoidal forcing. 
Furthermore, the absorbed energy averaged over the pe¬ 
riod of the perturbing field or over the orbital frequency 


of the electron is approximately similar to a resonance 
curve with one dominant frequency with finite peak at 
to = to a ; other small resonance peaks occur to the left or 
to the right of the major resonant frequency. 

The Lorentz dispersion model of the light-atom inter¬ 
action assumes that the electron is subject to Hooke’s 
force and the force due to the oscillating electric field 
of the light. Interestingly, even if our initial assumption 
is an electron in circular orbit around the nucleus, we 
still obtained the same forced harmonic oscillator equa¬ 
tion as that of the standard model. We also obtained 
the same resonant frequency, though the actual peak is 
at a frequency slightly smaller than ujq. But what is new 
is that even though we did not put a damping term in 
our harmonic oscillator equation, we still obtained a fi¬ 
nite energy absorption at the resonant frequency w. We 
also computed the electron’s angular momentum vector 
and showed that its tip traces rosette patterns similar to 
epicycles [25], 

In the future work, we shall extend our work to the in¬ 
teraction of the hydrogen atom with elliptically polarized 
radiation. 
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